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Abstract. We analyze the structure of the Fundamental Measure Theory for the 
free energy density functional of hard sphere mixtures. A comparative study of the 
different versions of the theory, and other density functional approaches, is done in 
terms of their generic form for the three-points direct correlation function, which 
shows clearly the main advantages and problems of the different approximations. 
A recently developed version for the monocomponent case is extended to mixtures 
of hard spheres with different radii, and a new prescription is presented to obtain 
the exact dimensional crossover of those mixtures in the one-dimensional (ID) 
limit. Numerical results for planar wall-fluid interfaces and for the ID fluid are 
presented. 
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1. Introduction 

In the recent developments of density functional (DF) approximations for the free 
energy of hard molecules the Fundamental Measure Theory (FMT), pioneered by 
Rosenfeld Q and explored in different directions by other authors [g, |j, U [| , stands 
at a prominent position both because of its peculiar functional structure and because 
of the combination of success and pitfalls obtained with its different versions. In this 
paper we analyze the generic mathematical structure of the different density functional 
approximations and discuss its implication for the generalization of a recent version 
of FMT to mixtures of hard spheres with different radii. 

The theories developed to approximate the interaction part of the free energy 
density functional <&[p] = P{F[p] — Fjd[p]), with (3 — (ksT)" 1 , beyond the simplest 
local density and square gradient approximations, are all based on a similar scheme: 
the non-local dependence on the density distribution p(r) appears through 'averaged' 
or 'weighted' densities defined as convolutions 



with several (k — 0,1,..) weight functions, u>fc(r — r'), which depends only on the 
relative position of two points. The density functional <&[p] is then written as the 
volume integral of a function of p(r) and/or the weighted densities Pfe(r). Such 




(1.1) 
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structure of the usual approximations for contrasts with the generic exact density 
expansion of the free energy, 

*M = \J drip(ri) J dr 2( o(r 2 )/(ri 2 ) 

+ i j dr x p(n) j dr 2 p(r 2 ) / dr 3 p(r 3 )/(r 12 )/(r 23 )/(r 31 ) + 0(p 4 ), (1.2) 

in terms of the Mayer function /(r) = 1— e - ^ r ' for any pair potential energy <^(r) (we 
have denoted = r.j — r^). The first term in (1.2) has the simplest structure of the 
density functional approximations, with a basic weight function Wq(t) proportional 
to /(r), used to convolute p(r) at two different points. However, the second and 
followings terms in ( |1.2j ) cannot be reproduced within the generic forms used in the 
DF approximations. In the diagrammatic expansion used in the theory of liquids §, 
the free energy is given by the irreducible diagrams (with Mayer function links), which 
cannot be evaluated as a single volume integral of simpler factors, while the generic 
form of the DF approximations includes only reducible diagrams with u>fc(r) links. 

The second and third functional derivatives of (fT^) give the exact series 
expansions for the pair and triplet direct correlations functions respectively, which 
have the corresponding density expansions, 

- C ( 2 >(n,r 2 ) = * 2 f [p ) =/(r 12 ) + /(r 12 ) / dr 3 p(r 3 )/(r 23 )/(r 31 ) + 0(p 2 ), (1.3) 
dp{Ti)5p{v 2 ) J 

and 

- c (3)( ri ,r 2 ,r 3 ) = — f j [P \ r = /(r 12 )/(r 23 )/(r 31 ) + 0(p), (1.4) 
5p{v 1 )dp(Y2)dp{Y 2 ) 

again made of irreducible diagr ams , with two and three open (not integrated) points 
respectively. The f irst term in (1.3) is reproduced by mo st D F approximations, but 
the first term in (|l.4| ), and hence the second term in (1.3) for arbitrary density 
distributions, cannot be reproduced within those DF forms. This contradiction 
between the exact forms and the usual choices for the mathematical structure of the 
DF approximations is forced by the computational cost of evaluating irreducible terms. 
The inclusion of kernel functions with more than two centers in the convolutions of 
p(r) would make it difficult to use the DF approximation for practical purposes. 

The weighted density approximation (WDA) [m, ||] gives a good example of the 
compromise between accuracy and practical computability in DF approximations. 
Although the advanced versions of this DF approximation use a density dependent 
weight w(r,p), which may appear to be beyond the scope of (1.1), for the present 
discussion it is fully equivalent to take the zeroth and first order contributions in the 
density expansion w(r, p) = Wo(r) + u>\{r)p + • • •. Then, the approximation for 
up to third order in p(r) is 

$wda[p] = ^ /drip(ri) J dr 2 p(r 2 )/(r i2 ) 

+ dr lP (n) J dr 2 p(r 2 ) J dr 3 p(r 3 )/(r 12 )[/(r 23 ) + Wl (r 23 )] + 0(p 4 ), (1.5) 



where the first term reproduces directly the first term in ( |1.2j ), while the second term 
tries to mimic the irreducible kernels as reducible combinations of Mayer function 
links, /(r), with a new link function, iUi(r), defined to recover the direct correlation 
function of a bulk fluid, i.e. to recover the exact expansion (fL3) up to first order in the 
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uniform density p(r) = p . This requirement leads to a function wi(r) which exceeds 
the range of the Mayer function with an oscillatory tail structure. 

Thus, the WDA for the hard sphere (HS) fluid may be implemented [|| to recover 
the direct correlation of the bulk fluid, c(rx2,pa) — c^ 2 )(ri, r 2 )| p ( r ) =Po , given by the 
Percus-Yevick approximation |(|, with the range of a HS diameter, c; however the 
function & 2 '{t\, T2), evaluated for a non-uniform density distribution, would exceed 
that range because of the convolution of the Mayer function with the oscillatory 
function ui\ (r) . Similarly, the zeroth order term of the triplet direct correlation 
t ffi ( ri , ra i ^3) would exceed the range = |rj — < a which the exact form 



( |l.4[ ) has for any combination of its variables. This being so, the WDA and similar 



approximations give a blurred representation of the exact non-local depende nce of $[p] 



on p(r), where the sharp step- like dependence on the HS Mayer function in (1.2)— (1.4) 
is replaced by an easier to compute but smoother functional. The extension of this 
type of DF to mixtures of HS with different radii R4 (i = 1, . . . , m) is difficult, because 
they require m(m + l)/2 weight functions to reproduce the Mayer functions fij(r) for 
any pair of molecules, and the 'extended' weights required to mimic the irreducible 
diagrams for any triplet (i,j,k = l,...,m) proliferate into an embarrassing large 
number of possible combinations for the functional 

The FMT also approximates for hard sphere systems in terms of convolutions 



(1.1) of p(r), but with weight functions having the range of the molecular radius, R, 
rather than the range of the Mayer function, a = 2R. The first advantage of such 
a choice is a natural extension to HS mixtures, since the number of different weight 
functions is going to be proportional to the number of species, and (at least in the 
original version) they only appear in a fixed number of linear combinations. The first 
one is the local packing fraction, which for a mixture of (i = 1, . . . , m) species of hard 
spheres with radii Ri and density distributions pt(r), is defined as 
m « 

r?(r) = / dr' Pi(r')e(^ - |r - r'|), (1.6) 
»=i •* 

with the Heaviside step function Q(x). This definition of ?y(r) is the direct extension 
to inhomogeneous systems of the usual packing fraction 77 = (47r/3) i?f p^, which 
plays a most relevant role in the best approximations for the equation of state of bulk 
fluid HS mixtures. Also, ?y(r) is a key ingredient in the exact density functional for 
hard rods pQ ] in one dimension (ID), and it has the appealing interpretation of being 
the probability that point r happens to be inside a sphere. Despite these facts, the 
local packing fraction cannot appear in any DF based on the Mayer function and its 



convolutions, like the WDA in (1.5), because all the weight functions have at least 



the range of two added molecular radii. The problem to include rj(r) in the functional 



depe nden ce of <&[p] is to recover the Mayer function in the lowest order terms in (1.2) 



and (1.3), from kernels with half the range. In the exact density functional for hard 
rods in ID, the function fij(x) = Q(Ri + Rj — \x\) is recovered as a convolution of the 
'molecular volume' Q(Ri — \x\) with delta functions, 

= \j Ax ' t e (^ - \ x ~ - M) + ®( r j -\ x ~ x '\) s ( R i - WD] > (i-7) 

which points to the use of normalized spherical delta function shells as 'molecular 
surface' weight functions to complement rj(r). 

In 3D the minimal set of weight functions with range Ri and Rj to recover the 



Mayer function spherical step at |r| = R4 + Rj, as in (1.7) for the ID case, requires a 
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scalar and a vector weight functions to define two sets of weighted densities: 

1 



rii(r) = 



dr , S(\r'\-R i )p i (r + v'), 



and 



Vi(r) 



4ttR- 



J dr'8(\r'\-R i )p i (r + r')^-, 



(1.8) 



(1.9) 



for each molecular species (i = 1, . . . , m) present in the system with density Pi(r). The 
structure of the density functional $[/?] at quadratic order in rj (r), rii(r) and Vj(r) is 
uniquely determined to recover the leading terms in (1.2) and (1.5), 



m = J < 



dr ( r?(r) ^ m(r) + 2vr ^ i?,i? 3 (i? 4 + i?,)[n 4 (r)n,(r) - v 4 (r) • v,(r)] 



+0(p 3 ), (1.10) 

The original proposal for the FMT |Q is consistent with ( 1.10 ) and, with 
the guideline of the scaled particle theory, recovers the full Percus-Yevick (PY) 
approximation for the direct correlation function of a bulk HS fluid. According with 
an general rule emanating from the FMT structure, the free energy density functional 
has as many additive terms as the space dimension, 



D 



For 3D hard spheres the three terms proposed by Rosenfeld were 



dr log[l — r](r)] 



i=l 



$ 2 [p] = 2ir RiRj{Ri +Rj) J dr 



rii (r)nj(r) - v,(r) 



v i( r ) 



*3,o[p] 



87T 2 ^ flf^^ 



drj n,-(r) 



1 - n(r) 
|n 3 -(r)7i fc (r) - 



/(r)-v fc (r) 



[1-J?(r)]2 



(1.11) 



(1.12) 



(1.13) 



(1.14) 



All the complexity ari sing from the mixture of m different HS species is reduced to the 
evaluation of r](r) in (L6) and three moments jW(r) = RfXi(r), with s = 0, 1, 2, 
of the averaged densities X{ — rx,,Vj. Thus, the theory may be applied even to 
polydisperse systems with a continuous distribution of the molecular radius [p"l|| . 

For one dimensional (ID) hard rods the ex act f ree energy DF, as found by Percus 
[|10|, is recovered from just the first term in 



1.12), 



with the obvious translation of 
r](r) and n^(r) to ID. For hard disks in two di mens ions (2D), Rosenfeld [[l2| proposed 
two terms with structures similar to ( 1.12| )— ( 1.13 ). The results obtained with this 
approximation for a monocomponent HS system were very good (better than those 
obtained with any previous DF approximation) for problems like determining the 
profiles of a HS fluid against a hard wall, thus showing that the non-local dependence 
obtained with the geometric measures (1.6), fll.Sj ) and (l.S) is a better representation 
of the sharp dependence of the Mayer function in the irreducible kernels of (1.2). 
However, the study of the HS crystal as a density distribution made of narrow 
peaks at the positions of the crystal lattice, which had been successfully achieved 
with the WDA J^, || and other DF approximations |l3| , gave fully unphysical results 
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with a negatively diverging free energy in the limit of infinitely narrow peaks. The 



problem was related to the overlap of three delta- function shells in (1.14), giving rise 
to integrable divergences. These singularities are harmless in the evaluation of the 
bulk fluid free energy and direct correlation function, but they notoriously show up 
when the free energy is evaluated for delta-function density distributions. In the case 
of a monocomponent HS system, the pathology of the original FMT was analyzed in 
terms of the dimensional crossover of the DF from 3D hard spheres to ID hard rods, 
as well as the zero dimension (0D) limit, i.e. narrow cavities which cannot contain 
more than one molecule [[l4j A new version of <J>3[p] was proposed to solve the 
problem |1^| , with the new ingredient of a 2-rank tensor- weighted density %, with 
cartesian components 

?f ,/3) M = 4^2/ dr '*(l r 'l + (1-15) 

a, (3 = x,y, z and where the index i is included to provide the obvious generalization 
to HS mixtures. This tensor-weighted density allows to write ^[p] in such a way that 
it vanishes for any ID distribution of HS, i.e. for any p(r) = S(x)S(y)pi(z), which 
should be fully equivalent to a ID system of hard rods with density distribution p±(z). 
The combined forms of $>i\p ] and ^[p] give the exact DF form in this limit, but the 
original term $3^] in ( |l.l4| ) spoils the agreement. The new version of the FMT for 
the monocomponent HS system reproduces the exact ID limit and gives an excellent 
description of the HS crystal, solving many of the qualitative problems of other DF 
approximations Jig] , 

The extension of this new version of FMT to HS mixtures amounts to taking 

* 8 M = Y^Rl^Rl J dr Ir ^gp, (1-16) 

with 

Vufc(r) = Vj • Tj ■ v k - njVi ■ v fc - Tr[TiTjT k ] + rij Tr[7^T fc ], (1.17) 

in terms of the tensorial contractions and trace of Vj(r) and Tl{r), for each component. 

However, such a direct extension to mixtures fails to reproduce the exact ID limit. 
The reason for this failure, the modification in the DF structure required to achieve 
that limit and the effect of such modifications in typical problems are to be analyzed 
here. In the next section, the origin and effects of the spurious divergences in the FMT 
density functional are analyzed in terms of the triplet direct correlation function and 
a new FMT density functional for HS mixtures is proposed to recover the exact ID 
limit. This new version, however, patches out this effect at the expense of creating 
new weaker but potentially harmful singularities which preclude using this functional 
for a free minimization. In section 3 we compare the results of this new version with 



those obtained with (1.16) and (1.17), and discuss the pros and cons of using each. We 



conclude in section 4 that the functional (1.16), (1.17), despite its defects, provides 



the best balance between accuracy and computational simplicity, and that it seems 
impossible to go any further without modifying the basic structure of FMT. 



2. FMT for HS mixtures 



The triplet direct correlation functions obtained for HS in three dimensions with any 
FMT density functional has a peculiar structure, very different from that obtained 
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Figure 1. Boundaries in the plane (r23,r3i), for the given value of ri2 = 0.4. 
The dash-dotted lines give the triangle inequalities \t23 — r"3i| < J"i2 < 7*23 +^3li 
for the possible values of these variables. The solid thick curves give the regions of 
common overlap between three spheres placed at ri, r2 and r3 while the solid thin 
lines give the boundaries of the simultaneous pairwise overlap. The HS radii arc 
of equal size (Ri = Rj = R^ = 0.5) (a) and of very different size (Ri = Rj = 0.5 
and R k = 0.1) (b). 



with the WDA or similar approximations. The most relevant form appears in the 
contribut ion a t zero order in a density expansion, from the functional derivative of 
$3[p] in ( 1.11 ), which has the generic form 

+ 0(p), (2.1) 

where i,j,k = l,...,m. The radii of the HS placed at ri, r 2 and r3 are Ri, Rj 
and Rk respectively, and for each possible choice of the three indices, i,j,k, there 
is a generic function Kijk of the relative vectors rj = rj — r, I = 1,2,3, whose 
moduli are restricted to be equal to the respective HS radii by the three delta- 
function shells. Higher order terms in the density expansion of c^, as well as the 
contributions to the triplet direct correlation function from the terms and $2[p] 

in (1.12) and (1.13), have one or more delta- function shells substituted by spherical 
step functions (i.e. Q(Ri — etc.) coming from the functional derivative of 

?7(r). These contributions are always regular smooth functions, which cannot produce 
spurious divergences but cannot reproduce either the exact step-like dependence of 
(Hjk (r l! r 2 ,r 3 ) = -fij{ri 2 )fjk(r23)fu{r3i) + O( i o), in terms of the Mayer functions, 
fij(r) — <d(Ri + Rj — r), and the relative distances ri 2 = |ri — r 2 |. 

It is only from the overlap of the three delta- function shells in (2.1) that the c^X 
term in a FMT may be able to reproduce such a sharp step dependence. But the same 
overlap is also the origin of the spurious divergences which may invalidate the use of 
the FMT for density distributions with strong inhomogeneities. The crucial point to 
balance these two aspects is the choice of the kernel function Kijk(r[, r 2 , r 3 ), which 
is precisely the difference between different versions of the FMT. The common aspect 

/"ON 

of all versions is that the range of the triplet correlation function, c^-((ri, r 2 , r^), is 
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restricted to points r 1; r 2 and r 3 such that the three spheres of radii Ri, Rj and 
-Rjti respectively centered at those points, have a common overlap. This geometrical 
boundary is always inside the exact range imposed by the Mayer functions product 
(which extends as far as pairwise overlaps between the three spheres exist) but, 
particularly in the case of mixtures of HS with very different sizes, it is much more 
restrictive than the exact one. In figure [l] we represent the boundaries in the plane 
(?"23, ?"3i) for a given value of m, both for HS of equal size (a) and of very different 
size (b) . The total hyper-v olum e obtained by integration over r 2 and r3 in the region 
with non-zero values for ( |2.l|) may be obtained as a function of the HS radii as 
while the equivalent integral for the exact product of Mayer 



,(3) 

FMT — ± H ±, 'j x ''ki 

functions gives 



r (3) 

exact 



16vr 2 



E?3 r>3 i d3 r>3 . d3 r?3 

+3R 2 R 2 Rl (z 1 Ri + R i i R i+ R k , Rk + Ri 



Rk Ri R 



■3 



(2.2) 



(3) (3) 

The difference between r^ xact and _T™ T corresponds to the FMT 'lost-cases' described 
in the OD approach to the FMT p5fl , i.e. those configurations of three delta- function 
peaks which cannot give contribution to 3?3[p] in the FMT scheme, but which should 
give a contribution to the exact free energy excess. When the three HS radii are 
equal the FMT contribution covers about 74% of the exact hyper-volume, so that 
the role of the 'lost-cases' is relatively mild, and its effect on the overall accuracy 
of the DF approximation for usual applications could be reduced by an appropriate 
choice of the kernel function in (2.1), which could increase the contribution of the 
points inside the FMT boundary to overcome the lack of contribution of the missing 
configurations. However, the ratio between T^-p and r^ act goes to zero in a mixture 
of very asymmetric HS, whenever the size of one species goes to zero while keeping the 
size of the other two. This is shown in figure || and corresponds to the effect which was 
obvious for the particular case in figure 0(b): the contribution of $3 [p] associated to 
spheres of very different size shrinks to a very small part of the exact contribution to the 
triplet direct correlation function. If the kernel function in (2.1) is fixed to reproduce 
a given equation of state for the bulk fluid mixture, the contribution of order piPjPk 
coming from ^[p] would require an enormous artificial enhancement of the weight 
of those configuration within the FMT boundary, in order to compensate the lack of 
contribution from the missing configurations. The use of such an approximation for 
inhomogeneous density distributions could lead to quite an erroneous evaluation of the 
excess free energy, with spurious sensitivity to the density distribution of the small 
spheres. This could be of particular relevance in the DF study of depletion forces in 
colloidal particles for strongly inhomogeneous density distributions. 

Beside the problem of its reduced range, the mathematical form (|2.l|) leads to 
the most serious problem of the FMT: the presence of spurious divergences at the 
boundary of non-zero values, arising from the tangency of two or the three spherical 
shells. This is clear if the integral over r is transformed into integrals over the three 
moduli r[, rU and r 3 , which can be directly performed thanks to the three delta 
functions in ( |2.l| ) . The boundary of non-zero values is therefore given by the existence 
of a triple overlap between the three spheres; since at that boundary the vectors r' 1; 
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Figure 2. Ratio between = 4tt 4 Rf R\ and r^ act [equation (Q] for 

Rj = Ri~> Rk (solid line) and for Rj = R^ < Ri (dashed line). 



r' 2 and r 3 are geometrically fixed, integration leads to 
<5 3 $ 3 [p] _ 4R i R j R k K ijk (T , 1 ,x f 2 ,T , 3 ) 



6pi(ri)dpj{r 2 )Spk(r3) [Z(r 12 ,r 23 ,r 31 , Ri, Rj , Rk)} 1 / 2 
The function Z(ri 2 , r 23 , ?"3i, -Ri, Rj, Rk) is the polynomial 
Z = (j£r§ 3 + iJ?i^)(r? a - r 2 2 3 + r§ x ) + (J^rfi + iZfo 2 )(r. 



O(p). 



2 

+ (-Rfe r 12 + ^i^j)( r 31 ~ r 12 + r 23) ~ J S ' 23 _ JL j' 31 ~ 1L k' 12 — ' 12' 23' 31 

arising from the Jacobian of the transformation. This polynomial vanishes over the 
whole boundary. The divergences of the terms in (2.1) and (2.3) may only be avoided if 
the kernel also vanishes over the boundary. The version of the FMT based on the exact 
OD limit for narrow cavities of arbitrary shape |Te| ] , fulfills precisely such a condition 
with a function ifj 7 - fc (r^, r 2 , r 3 ) proportional to \Z(r 12 , r 23 , r 31 , R h Rj, Rk)] 1 / 2 . This 
removes the divergences in (2.3) while at the same time saves the step- like dependence 
between the inside and the outside of the FMT boundary. However, such a kernel 
function is not separable in powers of the cartesian components of the vectors rj , so 
that the free energy cannot be evaluated in terms of simple weighted densities, like 
(1.1), making it very costly to use it for practical applications. 

The original FMT version (1.14) corresponds to building a kernel function on the 
basis of the bulk fluid direct correlation (from the integration of c\^ k with respect to 
one of the positions, as in (1.3)) and to get an easy computation of $[/?] in terms of 
the scalar and vector weight densities ([h^) and (1.9). In that case the kernel function 
has to be decomposable in terms which are linear on each cartesian component of 
the vector variables. Within thes e co nstraints the unique choice was ( 1.14 ), which is 
equivalent to taking in (|]l]) and (2J3) the kernel 

l 1 J. I A 1 ■ 1. ' ' a. " I I I r. 



31 ' 

4„2 



'31 



Rfr. 



R 4 r 2 



n k r 12 



2 2 
' 1 o^o 



(2.3) 



(2.4) 
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However, such a choice for the kernel function does not eliminate any of the divergences 
at the boundaries of the tangent spherical shells. The strongest effect appears when 
Ri — Rj — Rk and ri = r 2 = r 3 , so that the three spherical shells overlap in their 
whole surface. In that case ( |2.5| ) takes the value — (127r) _1 , and there is a diverging 
negative contribution to the free energy. The divergence is integrable and it contributes 
to give the exact p 3 term in free energy (1.2) for uniform bulk fluid, but in a density 
distribution with narrow peaks (like in a OD cavity or in a crystal phase) the negative 
divergence goes to ^s[p] in (1.14) and the approximation becomes useless. Also, in the 
crossover from the 3D hard spheres to t he ID system, there is a spurious contribution 
of the divergence in $3[/j], as given in ( 1.14 ), while the exact ID density functional 
would be recovered wit h a v anishi ng ^s \p] for such density distributions. 

The FMT version duel ) and (|L17|) was developed to recover the exact ID limit 
for monocomponent HS systems It corresponds to choosing the kernel 



K l]k 



I 

167T 



1 



Ri R 1 



1 - 



Ri Rk 



1 - 



Rk Ri 



(r: 



x r, 



OR; Rj R k f 



(2.6) 



The presence of a quadratic dependence in the cartesian components of the vector 
variables implies the need of the tensor weighted density, with rank two, defined in 
(1.15), together with the scalar and vector weighted densities (f[[|) and ( jl.9| ) used in 
(1.14). The kernel has two terms with different geometrical structure: one goes with 
the triple vector product • (r' 2 x r 3 ) which becomes zero whenever the three vectors 
r[, r' 2 and r 3 become coplanar. This happens always at the boundaries of the triple 
overlap of the three spherical shells in (2.1) so that this contribution to the kernel 
alone would get rid of all the spurious divergences generated by the vanishing of the 
denominator in (2.3). Unfortunately, a kernel function based only on this term would 
not reproduce the correct behavior for c^ 2 \r,p) in a homogeneous fluid. The other 



contribution in ( |2.6| ) becomes zero whenever two of the three vectors r[ become equal, 
which suppresses the strongest divergences of c^ 3 '(ri, r 2 , r 3 ) produced when two delta 
functions shells, of equal radii, overlap their full surfaces. In the monocomponent case, 
this is enough to make the full contribution of "^[p] to vanish in the ID limit, because 
three spherical shells of equal radii and with the centers along a straight line, cannot 
have a triple overlap unless two of them share the same center. 

Thus, in a monocomponent HS system the FMT version with ( 1.16 ) is free of the 
strongest divergences in (ri, r 2 , r 3 ), produced when ri — r m and when the three 
points ri,r 2 ,r3 are along a straight line (so that the spherical shells have a common 
axis). However, there are still weak divergences at the boundary points between 
overlapping and non-overlapping shells. These divergences would appear in for 
very peculiar density distributions, with three delta function peaks, normalized to 
contain less than one molecule among the three, and separated by distances such that 
the surfaces of molecules centered at the three sites would intersect in a single point. 
Of course such a density distribution would be very unlikely to appear unless the 
external potential acting on the system has been tailored on purpose, and in that case 
the error made by this FMT version would be of the opposite sign to the error made 
by ( 1.14 ) in the OD limit. I nste ad of a negative divergence of the free energy (from the 
negative sign of the kernel (2.5) in such a limit), there would be a positive divergence of 
the free energy because the kernel (2.6) is always equal or larger than zero for coplanar 
vectors. A positive divergence in &3[p] for a very peculiar density distribution would 
produce little harm because the minimization process of the free energy in the density 
functional formalism would just avoid such distribution functions. The problem is 
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Figure 3. Three coaxial spheres of different radii overlapping along a 
circumference. This situation is impossible for spheres of the same radius, unless 
two of them overlap their whole surfaces. 



much more innocuous than a negative divergence of the free energy, which would 
invalidate the formal use of $[/?] to minimize with respect to any possible function 
P(r). 

In the case of HS mixtures with different radii, the kernel (2.6) is less efficient 
than in the monocomponent case in the regularization of the FMT. The exact ID 
limit would still require that $3[p] vanishes whenever all the points with non-zero 
density are along a straight line, but the kernel (2J3) does not vanish in the case of 
three coaxial spheres with a triple overlap along a circumference (as illustrated in 
figure ||). To recover the exact ID density functional from a mixture of 3D hard 
spheres confined along a straight line requires a kernel Kijk(r[, r' 2 , r 3 ) which vanishes 
whenever the three vectors (with respective moduli Ri, Rj and Rk) go from a common 
origin to three points along a straight line. The simplest geometrical construction with 
this requirement is 



ijk 



(r' x x r 2 + r 2 x r 3 + r' 3 x r[ 



WirRiRjRk 



'3 



Ri + Rj Ri 



Rj + Rk Rj 



Ri 



(ri x r 3 ) 



Rk + Ri Rk + Rj 



(rj x r' 2 ) 
Rk Rj 
(r 2 xri) 



Rj Ri 



(2.7) 



r' 2 ), 



(2.8) 



Ri Rk 

The way to obtain this kernel is to multiply the vector 

ri x r 2 + r 2 x r 3 + r 3 x ri = (r' 3 - r 2 ) x (ri 

which obviously vanish whenever the three vectors are coaxial, by some linear 
combination of the three cross products in such a way that the contribution to the 
triangle in the virial expansion of the direct correlation function is recovered. Such a 
thing cannot be accomplished if the coefficients are just functions of the radii, but it 
can if we allow for a dependence on the dot products of the vectors. 

The existence of cubic dependence on the vectors in ( [2.7| ) forces to use a new 
3-rank tensor measure, namely A4i, whose cartesian components are given by 

' r dT 'S(\r'\-R i )p i (r + r') r ^l. (2.9) 



Mf aA7) (r) 



AnRi 



R: 



In terms of this and the former measures, ^[p] will have the same shape as in (1.16), 
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but with <^jfc(r) in (1.17) replaced by ip^iv) + At^^r), where 



Aip 



ijk \ 



2RAE4 — Rk) 



-(vj • Mj -.Tk-Vi-% 



(2.10) 



Ri{Ri + Rj){Rj + Rk) 

the symbol ':' denoting the contraction of two indices. An important difference arises 
in this new version: the coefficients multiplying the measures are no longer a product 
of powers of the radii, but more general rational functions of them. This has the 
computational disadvantage of not allowing to express the functional in terms of a 
finite number of moments of the weighted densities, as in any previous FMT version. 
This is particularly important in applications to polydisperse system, where having 
moment dependent functionals dramatically simplifies the calculations [ jTT| . 

A different and more formal drawback of the new kernel is the fact that it cancels 
the most singular terms in the former kernel (those that preclude it to recover the 
one-dimensional functional for aligned configurations of the HS) at the expense of 
changing the sign of the weak divergences for some non-aligned configurations, at the 
boundary between overlapping and non-overlapping shells. These divergences were 
always of positive sign for the kernel (2.6), but for some configurations they become 
negative with the kernel (2.7). This can be explicitly seen if we write the kernel (2.7) 
as 



K^ k =[(if 3 -x( i )x(y i -xi l )].V ijh (x{,i 



2- 



(2.11) 



and now rotate r' 2 slightly away of coaxiality while maintaining the three vectors in 
the same plane and with common origin (see figure |^) . If the displacement is given by 
the vector e (|e| <C Rj) then R 3 = \r' 2 + e| = Rj + 2r' 2 ■ e ■ 



Of 



0(|e| 2 ), i.e. 

(2.12) 

for coaxial vectors) 



Replacing r' 2 by r' 2 + e in (2.11) yields (remember that 

K ijk = [(v' 3 - ri) x e] • V l]k (r[,r' 2 ,r' 3 ) + 0(|e| 2 ). (2.13) 

Clearly, if e fulfills ( 2.12| ) and the r.h.s. of ( [2.13 ) does not vanish (which happens, for 
instance, if r' 2 and e are coplanar), so does — e, so the kernel may have both signs, as 
stated. 

A negative singularity is formally the worst defect a free energy density functional 
can have, because it implies that the absolute minimum of the free-energy functional 
is minus infinity, and it is reached for very singular density distributions. In practice, 
however, this defect can be overcome with a restricted minimization of the free energy 
within a functional family which does not include those singula r dist ri butio ns. Thus, 
the strong negative divergence of the original FMT proposal ( 1.11 )-( 1.14 ) for delta 
function density distributions, which precluded its use for crystals and narrow 0D 
cavities, does not interfere with its use (and very good results) for systems with 
planar symmetry, when the free energy is minimized with respect to planar density 
profiles pi(r) — pj(z) , ind ependent of x and y. The remaining neg ative divergences 
in our new version ( |2.13| ) are by far much weaker than those in (1.14), and they 
would be harmless for almost any virtual application of density functional theory with 
the variational minimization of the free energy constrained to families of functions 
(uniform in one or two variables, like in the adsorption at walls or in pores; periodic, 
like in freezing; etc.). This notwithstanding, we will show in the next section that 
the actual quantitative difference between this new functional and the former one 
(equations (1.16) and (1.17)) in applications to standard problems is negligible, the 
former functional being far simpler and easier to apply to polydisperse mixtures. 
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Figure 4. Maximum, with respect to the radii ratio R[/R s , and the relative 
packing fraction of the small segments, x s = p« /w - of th e rela tive error A^/^id of 
the one-dimensional reduction of functional (1.16) and (1.17) for a binary mixture, 
versus the total packing fraction, rj, of the fluid. The values of the pair (x s , Ri/R a ) 
are given for three points of the curve. It can be seen that the maximum error 
occurs for x s m 0.2 and Ri/R s m 4 for most values of rj. 



3. Results 



As stated in the previous section, the functional ( 1.16j ) and (1.17) does not recover the 



exact free energy of a uniform hard rod mixture when it is reduced to a one-dimensional 
density profile. The expression of the reduced functional is rather cumbersome, but 
in order to estimate this deviation we will simply evaluate the difference with respect 
to the exact value for a binary mixture. Let us denote respectively $ and ^id the 
reduced free energy density (in units of ksT) and the exact free energy density, and 
define A$ = $id — The relative error A$/$id has, for any packing fraction 77, an 
absolute maximum as a function of the relative packing fraction of the small segments, 
x s = rj s /ri and the radii ratio, Rl/R s - This maximum value is plotted vs. rj in figure^. 
In figure |^ we also plot the relative error A$/$id corresponding to rj = 0.5 both vs. 
x s at fixed Rl/R s (a) and vs. Rl/R s at fixed x s (b), for the values corresponding 
to the maximum relative error. Notice that this error remains smaller than 2% up to 
rj = 0.8, and it is never larger than 10% even for packing fractions as high as 77 = 0.95. 
The conclusion we extract from these figures is that, despite having an incorrect 



dimensional crossover to ID, the functional (1.16), (1.17) produces results of high 



accuracy in this limit. A second test can be obtained by applying this functional and 



the one with the correction (2.10) to obtain the density profiles of a binary mixture of 
hard spheres both, near a hard wall and within a slit. The parameters of the system 
have been chosen in order to obtain the most prominent differences. The results are 
shown in figure [| for the case of the hard wall, and figure fj] for the case of the slit. 
We see again that differences are negligible. 



In view of these numerical tests we conclude that functional (1.16), (1.17) is 
preferable to the corrected one, despite its inability to recover the one-dimensional 
free energy exactly, simply because its structure is far simpler and more suitable to 
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Figure 5. Rela tive error A<I>/<I>id of the one-dimensional reduction of functional 
( L. 16 ) and (1.17) for a binary mixture at total packing fraction rj = 0.5, versus 
the relative packing fraction of the small segments, x s = r) s /ri, for a radii ratio 
Rl /R s = 4.22 (a) as well as versus the radii ratio, Ri /R s , for a value of x s = 0.211 
(b). 
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Figure 6. Density profiles of a binary mixture of hard spheres of radii f? „ = 1 
and Ri = 5 near a hard wall, as obtained both with functional flLlq ), (1.17) 
only (dashed line) and with the correction (2.1C) (solid line). The densities are 



multiplied by the volume of each type of particle. The bulk packing fractions of 
the small and large spheres are, respectively, r) s = 0.2 and rji = 0.3. The difference 
between the results of both functionals can hardly be seen in the plotting scale. 



be applied to polydisperse mixtures (see for instance |11|). 
4. Conclusions 

The main conclusion of the analysis we have performed on the construction of a 
fundamental measure functional for mixtures of hard spheres is that, as stated in the 
title, we appear to be 'close to the edge' of fundamental measure theory, in the sense 
that its internal structure seems almost exhausted. There are two main drawbacks 
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Figure 7. The same as in figure |6| but for a slit of length 20-R s . The system 
is in chemical equilibrium with a reservoir with packing fractions rj s = 0.2 and 
r}i = 0.2. Although the difference is now more visible, it is still negligible. 



that seem unavoidable within the present scheme, both related to the structure of 
the third term $a[p] in ( 1.11 ). The first problem is common to all the FMT versions 
and it concerns the existence of 'lost cases', i.e. the fact that configurations of three 
spheres with pairwise overlap but no triple overlap do not make any contribution to 
the functional. This reflects in that the lowest order of C3 is incorrectly predicted (it 
vanishes for the lost cases, where it should still be 1). The obvious way to repair 
this defect is to introduce two point measures joining two 'halves' of Mayer functions. 
This would be a qualitative departure from the FMT structure, which would recover 
the exact C3 at low density, but the increase in computational complexity would 
be considerable. Another way to circumvent this problem, while keeping the FMT 
structure and low computational cost, would be to replace the delta functions in 
the definition of $3[p] by some other functions having a 'tail' extending beyond the 
radii of the particles (and presumably vanishing when particles no longer overlap). 
Although we cannot definitely exclude this possibility, we have tried several functional 
forms without success. This seems to be a too drastic change within the fundamental 
measure scheme, so much that adding such tails would 'break down' the functional 
at some other point. In particular, the dimensional crossover (to the exact ID limit 
and to OD cavities, strongly related to the density distribution in crystals) is a very 
stringent, hence fragile, requirement for an approximate functional, so virtually any 
modification can spoil it. 

The existence of these lost cases forces to overweight the other configurations in 
order to compensate for the nonexistent ones and guarantee a reasonable equation of 
state. But as figure [l] and figure || show, this compensation introduces a dramatic 
bias in the case of mixtures of very dissimilar spheres. This may, at least in principle, 
have important consequences for the depletion effect in these systems and thus may 
affect phase behaviour, at least in its use to describe crystalline phases. Although in 
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the application to planar density profiles, the effect would be strongly reduced by the 
averaging over the transverse directions. 

The problem of the lost cases seems to be inherent to the fundamental measure 
structure with a notable exception: parallel hard parallelepipeds Q|. The peculiar 
shape of these bodies makes pairwise overlap and triple overlap to be equivalent 
conditions, so there are no lost cases for this particular form of the particles, and 
hence the lowest order in C3 is exactly recovered. 

The second problem, a lso associated to $3 [p] , depends on the choice of the kernel 
function Kijk(r[, r' 2 , r' 3 ) in (2.1), which marks the difference between the FMT versions 
proposed in the literature. The balance here comes between the computational cost 
and the existence of spurious divergences (of positive or negative sign) in the free 
energy. The original form proposed by Rosenfeld factorizes in terms of the scalar and 
vector weighted densities, but the negative divergences show up in the reduction to 
ID and to OD distributions. Simpler versions (for the monocomponent case) using no 
kernel function and only the scalar weighted density || give much stronger divergences 
(showing up even at the bulk fluid direct correlation function c(r, po)) but always 
with positive sign, so that in the minimization of the free energy those configurations 
contributing to the divergences are avoided, e.g. creating spurious kinks in the density 
profiles of a wall-fluid interface, but without the qualitative breakdown of (2.5) for 
strongly confined systems. 

In contrast, the FMT version based on the exact reduction to the OD limit |Q 
would have no divergences of the free energy at all, but it has a non-separable kernel 
with enormous computational cost (implying the multiple integrals over the densities 
at three different points). The proposal made in jlj| and used by Groh and Mulder 
p7[ of taking ~ [r^ • (r' 2 xr' 3 )] 2 would also kill all the divergences at the boundary 
between overlapping and non-overlapping spheres, and be separable in terms of the 
scalar, vector and a rank-two tensor weighted densities. However, this form gives a 
too smooth boundary behaviour for (2.1), since the kernel vanishes as the function 
Z in the square root of denominator in (^), while the step- like dependence would 
require a non-separable kernel proportional to Kijk ~ • (r' 2 x rg)|. The alternative 
kernel ( |2 .6| ) and DF (1.16), proposed and successfully used ]l6| for monocomponent HS 
systems, achieved the separability (in term of the same set of weighted densities) and 
the sharpness (with the exact contribution to the triangle diagram for the bulk fluid) , 
at the price of unleashing some weak divergences of positive sign of the free energy for 
some peculiar configurations of non-aligned molecules. The effect of those divergences 
would be extremely difficult to observe and the quality of that DF approximation 
seems to be limited mainly by the quality of the Percus-Yevick equation of state for 
the bulk HS fluid, and the associated condition of having a direct correlation function 
with the range of the hard core diameter Jig] , which is also intrinsic to the FMT. 

The extension to HS mixtures imposes more severe restrictions, since the 
form (2.6) for HS of unequal size has (positive) divergences for peculiar aligned 
configurations which were not present for equal sized HS and which give stronger 
contributions than those of non-aligned configurations. Their main drawback is to 
spoil the exact reduction to the ID limit, which was one of the most remarkable 
achievements for the monocomponent case. We have introduced here a new kernel 
form (2.7) which would recover that dimensional reduction for HS mixtures, at the 
price of introducing a new rank-three tensorial weighted density and creating some 
weak negative divergences for some non-aligned configurations. Besides, the radii- 
dependent coefficients of this new kernel are not simply products of powers of the radii, 
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but more complicated rational functions, so it turns out to be computationally more 
involved in applications to, for instance, polydisperse mixtures. Other possibilities 
remain for a candidate kernel, like squaring the one we propose, conveniently adjusting 
the coefficients; however this forces the introduction of very high rank tensor weights 
and removes the discontinuity of C3, thus creating both computational and structural 
problems. 

Altogether, the systematic improvement of the free energy density functional 
within the FMT seems to be frustrated by incompatible requirements, deeply rooted 
in the structure of the theory. Thus, one has to choose a functional form on the basis 
of its performance for the required practical use. The original proposal by Rosenfeld is 
quite accurate for systems with planar symmetry, for which the imp rovement obtained 
with the newer versions is purely marginal. The FMT version ( 1.16] ), with a moderate 
i ncrea se of the computational cost due to th e presence of the tensorial weighted density 
( 1.15 ), would give similar results as ( 1.14 ) for planar systems but it is also good for 
crystalline phases or other strongly confined density distributions. In the ID limit 
it does not recover the exact DF for arbitrary distributions, but making a rather 
small error unless the ID system is very close to the close packing limit. The extra 
computational cost of the FMT version ( 2.1C ), with a rank-three tensor and non- 
additivity of the weighted densities in terms of a few moments of the radii, would only 
be worth for problems in which the exact ID reduction is crucial and the spurious 
negative divergences of this free energy DF are made innocuous by the parametrization 
of the density distribution. It seems very unlikely that these two constrains appear in 
a practical problem, so that probably (1.16) is the best practical choice. 

To conclude, we have to point that the problems and limitations of the FMT 
density functionals, which we have explored here, should be considered in the right 
perspective, as we are asking the FMT functionals to perform well in problems which 
would be out of question for any other type of DF approximation. The exact (or 
nearly exact) reduction from three dimensional density distributions to the ID limit, 
and the accurate behaviour in very narrow cavities, which opens the description of 
the crystalline phases with unconstrained density distributions, arc fully out of range 
for the WDA or other DF approximations. The presence of spurious divergences 
in the FMT is the price to be paid for the sharpness of its non-local dependence, 
compared with the blurred dependence of the other DF approaches. The taming of 
those divergences to make them innocuous has been the 'holy grail' of the workers 
in the fi eld, st artin g from Rosenfeld's choice for the combinations of n^(r) and v^(r) 
in ( 1.13| ) and ( 1.14 ). The 'earthly' requirements of practical computability make that 
'grail' unreachable, but the new versions explored here seem to be close enough to it for 
most practical purposes, providing by far the best density functional approximations 
for systems of HS mixtures. 
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